Rough edges in quantum transport of Dirac particles 
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We consider Dirac particles confined to a thin strip, e.g., graphene nanoribbon, with rough edges. 
The confinement is implemented by a large mass in the Hamiltonian or by imposing boundary con- 
ditions directly on the graphene wave-functions. The scattering of a rough edge leads to a transverse 
channel-mixing and provides crucial limitation to the quantum transport in narrow ribbons. We 
solve the problem perturbatively and find the edge scattering contribution to the conductivity, which 
can be measured experimentally. The case of Schrodinger particles in a strip is also addressed, and 
the comparison between Schrodinger and Dirac transport is made. Anomalies associated with quasi- 
one dimensionality, such as Van Hove singularities and localization, are discussed. The violation of 
the Matthiessen rule is pointed out. 



Recent discovery of graphene [2, [2] has spurred much 
interest. The honeycomb structure of the graphene lat- 
tice implies that carbon orbitals are sp 2 hybridized, leav- 
ing one free (p z ) electron orbital per atom. Unlike the 
majority of systems in condensed matter physics, at low 
energies these electrons effectively obey Dirac equations 
of motion for massless fermions |3j , with the Fermi veloc- 
ity being the speed of light. Another example is a d-wave 
superconductor whose electrons behave relativistically in 
the vicinity of nodal points 0] • The "relativistic" char- 
acter of graphene implies that its properties should be 
essentially different when compared to more traditional 
systems, whose electrons obey Schrodinger equation @- 

The central problem addressed in this Letter concerns 
transport properties of narrow graphene ribbons (strips). 
In general, the transport in thin and narrow structures 
is limited due to the restrictions imposed by quantum 
mechanics [6|. Classically, a beam of particles may be 
collimated as finely as desired so that it does not inter- 
act with the boundaries of a sample. As a consequence 
- in absence of ordinary impurity scattering - the mean 
free path of an electron is infinite, leading to vanishing 
resistivity. Due to the quantum mechanical uncertainty, 
however, a beam of particles must have a finite width: 
as the lateral dimensions of a sample shrink toward and 
below this width, the transport exhibits a crossover from 
the quasiclassical to the purely quantum regime, gov- 
erned by boundary scattering [6|. This is true both for 
Schrodinger and Dirac particles. We expect our results to 
have broad implications for graphene based electronics. 

We use the Kubo formalism Q to derive the expression 
for the conductivity. Similar approach was used in Refs. 
B S EH II , IH j although these authors were interested 
in other issues. Alternatively, the conductivity can be 
derived from the dielectric function 1131 . or equivale ntly 
3 from the Landauer formula [H, HI, EE El El, Hf- 
The quantum limitations to the transport properties of 
gra phene strips with smooth edges were considered in Ref. 
[15j : obviously, this is a different problem, exploring bal- 
listic transport, in contrast to the present paper, where 
the electrons are subject to a weak random scattering 



(impurities, phonons, etc.), and have a finite mean free 
path. The transverse channel mixing due to rough edge 
scattering in graphene junctions was numerically exam- 
ined in Ref. [2fJ, however, this study too is centered on 
ballistic transport. 

The confinement of Dirac fermions is a subtle issue due 
to the Klein paradox 2l|. This problem can be circum- 



vented by using a large mass term outside the strip in lieu 
of a static potential [22J . We employ this general method 
and refer to strips with such boundaries as the Dirac 
strips or ribbons. The boundary conditions in graphene 
are sensitive to the details of the edge l3|, and can also 
be formulated on the microscopic level [23( . For definitc- 
ness, we consider the special case of a graphene strip with 
metallic armchair edges, i.e., armchair strips/ribbons, 
and assume the variations in its width do not change 
the metallicity at the boundary 0, [3] ■ The spectrum of 
metallic armchair strips is exceptional since it contains 
a band with zero transverse momentum. Consequently, 
the armchair strips are anomalously good conductors, es- 
pecially in the vicinity of the Dirac point, where all the 
other boundary conditions yield a finite gap. At a finite 
chemical potential these differences subside but the afore- 
mentioned zero-band remains special as it does not get 
scattered from the edges 25] . Our boundary conditions 
assume "clean" edges, i.e., the dangling carbon bonds at 
the edges are inert and do not absorb environmental im- 
purities that might cause the current redistribution j2?| 
or localization of charge carriers 28] . 

The scattering matrix due to rough edges is Q 



V cx [\(x){y,p y },Hy 



(1) 



where x and y are the propagating and the trans- 
verse direction, respectively, H y is the transverse part 
of the Dirac Hamiltonian (straight edge), and X(x) 
parametrizes the width profile. In contrast to the 
Schrodinger case, the matrix Eq. |T]) is not separable. 
The matrix elements for a single rough edge can be found 
in Ref. [29| , and it appears that the only way to proceed 
requires inverting dense matrices of n c x n c size, with n c 
being the number of transverse channels. To overcome 
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FIG. 1: A thin strip of width d, whose bottom edge is rough 
(the thick lines represent the physical edges). A buffer layer 
(light gray) with large, but finite mass is attached to the rough 
edge of the sample. 

this sensitivity of the Dirac fermions to rough edges, we 
implement varying width boundary through an alterna- 
tive construction, which mimics the confinement of the 
Dirac fermions. In place of an infinite mass jump at the 
edge, we allow the mass parameter to change in two steps. 
At the edge the mass becomes large but finite, followed 
by the infinite mass jump (Fig. [J). This buffer layer of 
variable width suppresses the wave-function outside the 
strip while it simultaneously allows us to recast the equiv- 
alent of the scattering matrix Eq. {T]) in a quasi- separable 
form, amenable to analytic treatment. 

The Letter is organized as follows: first, we consider 
Dirac states of a thin strip, and derive the Kubo formula 
for these Dirac states. The Kubo formula is tested on 
a case of a smooth finite-width strip, whose resistivity is 
needed in what follows. Next, the buffer layer is intro- 
duced and the corresponding matrix elements that mix 
transverse channels are found. The exact mass parame- 
ter appearing in the matrix elements is determined by de- 
manding that the first-order perturbative corrections to 
the conductivity from a buffer layer with a straight edge 
are in agreement with the physical result. The effects of 
rough edges appear in the second-order perturbative cor- 
rections to the energy. This leads to our main result, the 
conductivity due to rough edge scattering. Finally, we 
find the conductivity of a thin strip whose electrons obey 
the Schrodinger equation The comparison is made 
between the conductivities of two identical strips, one 
containing Dirac, and the other Schrodinger particles. 

The 4-component Dirac equation 

(v 7 o7 • P - M)* = e*, ( 2 ) 

with 70 = diag(er 3 , cr 3 ), 71 = diag(-MT2, icr 2 ), and 72 = 
diag(i(7i, ici), has solutions in form of plane waves 



ttmC,a(r) = e 



ikx 



„(-) p -iq m y 



(3) 



In our notation, k is the momentum in the propagat- 
ing direction, and q m is the transverse momentum. The 
boundary conditions discretize the latter, yielding q m — 
(in + l/2)n/d, and q m = m-n/d in the case of a Dirac, 
and an armchair strip respectively, m is taken nonnega- 
tive. £ is ±1 corresponding to conductance and valence 
Dirac cones. The elements of spinors ut™ £ a are deter- 
mined through the boundary conditions [3[ and the Dirac 



equation Eq. ([2]). Index a = 1,2 denotes two decoupled 
valleys for a Dirac strip, and two orthogonal states which 
admix the valleys for an armchair strip. The energy of 
state Eq. ^ is £fe m j = —u + ^v^Jk 2 + q 2 ^, with fj, being 
the chemical potential (equivalent to the gate voltage). 
The unperturbed propagator is given by G™^ b a (iuJn', k) = 

Gmla( iw n] k)S mn S^S ab = 5 mn 5^5 ab / (lLO n - € kmi ). 

The (longitudinal) current operator for a Dirac particle 
is j x {^) = evif> (ry 7071^ (r) . Starting with the Kubo 
formula, and invoking the Lchmann representation for 
the propagtors 0, we find the one-loop contribution to 
the dc conductivity for Dirac particles 



Ode 



2nd 
(ka 3 



Y, ^.(o,*K$(o,*)x 

mQ^ab nrjpcd 



(4) 



-2 Im G™^ b a (E + iS; k) is the spectral 



Reve,A^ b a (E,k) 
function matrix. 

The Kubo formula Eq. |QJ reproduces the Dirac point 
conductivity, both the Gaussian [j| and the universal val- 
ues [3(|, depending on the spectral function used. We use 
the Eq. Q to find the conductivity of a straight strip of 
width d, where particles weakly scatter with an average 
lifetime r = l/(2r). A straightforward calculation yields 
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2tt dT 



Y 2 



(5) 



where Y m = vq m /\fj,\, and the star indicates summation 
over all m-s with Y m < 1. When a strip is wide (or /1 
large), the conductivity converges to a — §-(fj,/T) re- 
gardless the edge type. 

We now implement the confinement of the Dirac par- 
ticles by means of a finite mass buffer layer. The ba- 
sic idea is that in a lattice system, it suffices for the 
wave-function to fall-off faster than the lattice spacing 
to effectively prohibit hopping between the sites at the 
edge and their (non-existing) neighbours. Hence, mass 
M has to be large, and its precise value will be deter- 
mined shortly. Since the buffer layer serves as an ex- 
tension of the strip exterior, this mass should be of the 
same type as the infinite mass confining the fermions 
to the strip 3l|. For a Dirac strip that is the chiral 
mass, and a buffer layer of variable width a(x) intro- 
duces a perturbation V — Mv 2 j , for < y < a(x). The 
armchair boundary conditions are instead recreated via a 
( -ia 2 

turbation for this edge is V = Mv 2 ^^. Using t = —172, 
the matrix elements to the leading order in a are 



'mass term" ; the buffer layer per- 



V, 



InQa 



km£ b 



[(fL< a rW(hrn (b T] , (6) 
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where summation over fx, v is assumed, and [a(x) 2 ] 
is the Fourier transform of a(x) 2 . Denning z km ^ = 
£(k + iq m )/\/k 2 + q 2 n , we have f£ n(1 = (-iq m (l + 
Zkm£)A ~ z km £,Q,0) T , and = (0,0, -iq m (l - 

z femj)> (1 + z kms)) T f° r Dirac strips, while for armchair 
strips f£ m(1 = (-iq m Zkm£,l,iQm, z m k(,) T \ and fkm& = 
(fkm£i)*- Importantly, the channel mixing matrix Eq. 
(J6]) is quasi-separable. 

Now, consider a buffer layer of a uniform and small 
width ao. Both M and ao enter Eq. ([6|) as free param- 
eters. It is, however, clear that the first order correc- 
tions to the conductivity will be proportional to Ma\, 
while, according to Eq. (O, they should be proportional 
to ao- The conclusion is that the product Mao has to 
be a well defined constant. This is rather natural: as the 
buffer layer shrinks, the suppression of the wave function, 
driven by the inverse mass, is enhanced. To find this con- 
stant, we recall that the change of the conductivity by an 
innnitcsimally thin buffer layer should satisfy: 



So- = a(d) — a(d — a) w a d<Jo(d)/dd 



(7) 



The right hand side of Eq. Q follows from Eq. © . The 
left hand side can be found from Dyson equation 

Cl _ vknCa G (0)\-lykspc G {0) fo\ 

y nCa\ 1 v ksnc y spc) v km£b y m£b' \° I 



where argument (iu> n ; k) is assumed. The self-energy fol- 
lowing from Eq. (JSj) is real, yielding the "effective" dis- 
persion 4 m? = e hmi +£Mv 2 a 2 ,q 2 n / (d^/k 2 + q?~) for both 
Dirac and armchair strips. The corrections to the con- 
ductivity caused by this dispersion are now substituted in 
(J7J) , implying that the mass in ([6]) has the channel depen- 
dent form M = l/(vaoY m Y n ). Thus, an infinitesimally 
thin buffer layer perturbatively - at the lowest order - 
reproduces the physical behavior expected from a layer 
of infinite mass. A metallic armchair ribbon needs spe- 
cial attention: the coupling for scattering into its lowest 
energy channel (Yq = 0) appears infinite, thus rendering 
the perturbative approach invalid. This is reflective of an 
anomalous character of such a state, stemming from spe- 
cial boundary conditions. Since this state is impervious 
to the mass term, we ask how it is affected by the change 
of the strip width exactly, i.e., we directly examine the 
matrix elements Vo m , Eq. ([1]). They all vanish and here- 
inafter we ignore the scattering to/from this channel. 

We are now ready to find the second-order self-energy. 
Due to the quasi-separability of the scattering potential 
Eq. (O, we use the following ansatz W n ^ a ,mib{iu n \ k) = 
(fLcaT^p-( iu; n;k)(f kmib Y. The self energy and the 



propagator equations are cast as 
,Mv 2 a 



W(ico n ;k) = (- 



2d 



) / dl (w fc _iW;_fe> X 



(9) 



tgW(iu} n ; l)[l - W{ioj n )g^{ioj n - l)]~H, 

Sl% = G [ %Jn m 5 a tab+ (10) 

G^UaW Wg^W(f mib )G ( X 
The projected propagator is defined as 

rn£a 

In ([5]), w(x) is the deviation from the average buffer 
thickness ao- For simplicity, we assume the white noise 
edge profile characterized by (wkW-k) = aH 2 , where a is 
the lattice spacing, and H 2 measures the rms fluctuations 
in the strip width. Combining the self-energy from Eqs. 
(0 [TO]) and Kubo formula ([J|, yields the leading order 
conductivity expressed universally as 



(12) 




This is the main result of this Letter. 

Unfortunately, the summation in Eq. Q12p is difficult 
to carry out analytically and must instead be performed 
numerically for each \x and d. An example is plotted in 
Fig. [2^- Note sharp drops (to zero) of the conductivity 
each time a new channel is open. With [i at the bottom 
of a newly open band, the density of states experiences 
Van Hove singularity (absent in higher dimensions). As 
the number of states available for scattering increases, 
so does the self-energy, and the "effective" lifetime of 
excitations goes to zero. In Eq. (fT2|) . this corresponds to 
the situation when the highest Y m is close to one. 

For a wide strip, or large [i, one can estimate the 
asymptotic behavior of the maxima in the conductivity 
Eq. (fT2"|) , i.e., its values just before the opening of a new 
channel. For a Dirac strip, the maxima converge to 



@max,D 



e 2 v^'H 1 ^ 



2n 



1/4 



(13) 



In the case of an armchair strip with the same width 
the result is 3 3/4 w 2.28 times larger. Although Fig. [2ja, 
shows approximately equal conductivity maxima at low 
chemical potential, as /i increases, so does the armchair 
strip become better conductor as compared to the Dirac 
strip. The conductivity converges slowly to Eq. (fl3|) . as 
shown in the inset of Fig. [2^. 

For comparison, we consider the case of a strip with 
"Schrodinger" carriers. The Fermi energy is denoted 
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FIG. 2: The conductivities of a thin strip with a rough edge 
in units of e 2 /h. We set d = 1, and w = 1 (Dirac, a), or M = 1 
(Schrodinger, b). The roughness is d 3 /aH 2 = 10 4 . In all three 
cases, the Van Hove singularities appear; they are equidistant 
for Dirac particles, while their distance progressively grows 
in the Schrodinger case. The conductivity(ies) of the Dirac 
particles falls off considerably faster (fi^ 11 ^ 4 ) than that of the 
Schrodinger ones (/x -1 / 2 ). The inset shows the asymptotics: 
the dashed lines correspond to to Eqs. (113[) and (|15|l . The 
conductivity of an armchair strip (a; green) is on average 
greater by factor 2.28 than that of a Dirac strip (a; red). 

by n, the same as the chemical potential of the Dirac 
case. The dispersion in an unperturbed system is ekm = 
—fj, + (k 2 + q 2 n )/(2M), M being the electron mass here. 
The formalism of Ref. Q is readily adapted to a thin 
strip. Assuming the same parameters as before, the con- 
ductivity is 



2vr aH 2 M(i< 



Y 2 



V^Y 2 



(14) 



and is plotted in Fig. [2b- As in the previous case, 
Y m = q m /-\/2MJi. Van Hove singularities are appearing 
here too; however, in contrast to the Dirac case, the sin- 
gularities are not equidistant, and are instead distributed 
as squares of integers, due to the parabolic dispersion. 
The other qualitative difference is reflected in the asymp- 
totic trend for the conductivity maxima. For wide strips, 
the conductivity peaks follow the envelope function 



Ad 2 



& max,S 



2tt 3aH 2 y/2Mjl' 



(15) 



The conductivity has a sharp d 2 dependence (compare to 
d 3 for thin films @]). 

There is a notable qualitative difference between 
the asymptotic behaviors Eqs. (fT5|) and (fT3f . For a 
Schrodinger strip, the major contribution to the conduc- 
tivity comes from the lowest bands, as these states scatter 
the least from rough edges. Accordingly, the sum in Eq. 
(fT5")) is asymptotically proportional to the number of open 
channels, and the overall conductivity is proportional to 
^-1/2 rp^ e Tjj rac states in the lowest bands are, on the 
other hand, the most susceptible to the edge scattering, 
and their contributions are largely suppressed. Hence, 
the conductivity Eq. (fT3")) decreases rapidly (~ /i -11 / 4 ) 
as the chemical potential increases. The high impact 



of the edge scattering on the Dirac particles is also ev- 
ident in the weak d 1 ^ 4 width dependence in Eq. (|13p . 
These peculiar power laws are the consequence of the 
conductivity-per-channel function in Eq. (fT2|) . 

In quasi one-dimensional systems, one must be mind- 
ful of inevitable localization [32j . Its effects are discussed 
elsewhere (2(| . Generally, one anticipates a certain local- 
ization length L for electrons. The results presented here 
should be valid for weak disorder, i.e. whenever the strip 
is shorter than L. 

In thin films, the scattering from a rough surface yields 
a channel dependent mean free path, whereas the impu- 
rity scattering is channel independent, resulting in the vi- 
olation of the Matthiessen rule [33| • In thin strips the sit- 
uation is similar, the channel mean free path l' m from Eq. 
(I12p is not proportional to the impurity mean free path 
of a corresponding channel l® n — vwl — Y^j '{dTkp), Eq. 
(O, hence the Matthiessen rule is violated here, too. 

In this Letter we have developed a perturbative ap- 
proach to the problem of the edge scattering in Dirac 
strips. The particle confinement is implemented in the 
manner that allows us to solve the problem analyti- 
cally. The conductivity is found to the leading order 
in the edge roughness. We also analyze thin strips with 
Schrodinger particles and compare the results to those 
for Dirac strips. In both cases, the conductivity de- 
velops Van Hove singularities; however, different disper- 
sions mean that the singularities for Dirac particles are 
equidistant in the chemical potential, while those for the 
Schrodinger ones are not. The higher sensitivity of Dirac 
particles to the edge roughness suppresses the otherwise 
most conductive channels, and accordingly, the (average) 
conductivity fall off as a function of the chemical poten- 
tial is much faster as compared to the Schrodinger case. 

This work was supported in part by the NSF grant 
DMR-0531159. 
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